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ABSTRACT 

Wave-induced  vessel  motions  often  determine 
whether  launch  and  recovery  operations  can 
proceed  safely.  For  launch  and  recovery  of  a 
small  craft  from  a  larger  ship,  wave-induced 
motions  of  the  larger  ship  will  influence  dy¬ 
namic  loads  on  the  crane.  The  motions  of  the 
small  craft  will  be  a  major  determinant  of  the 
safety  of  onboard  personnel.  This  paper  exam¬ 
ines  wave-induced  motions  during  launch  and 
recovery,  including  the  effect  of  hydrodynamic 
interactions  between  the  two  vessels.  The  pres¬ 
ence  of  a  larger  ship  is  shown  to  have  a  major 
influence  on  motions  of  the  small  craft.  Sample 
polar  plots  demonstrate  the  variation  of  rele¬ 
vant  motion  parameters  with  speed  and  head¬ 
ing,  and  can  be  used  as  the  basis  for  operational 
guidance. 

INTRODUCTION 

Launch  and  recovery  of  small  boats  from  ships 
is  a  vital  part  of  naval  ship  operations.  Interest 
in  launch  recovery  among  navies  is  currently 
very  high  for  a  variety  of  reasons.  Interdic¬ 
tion  operations  have  become  increasingly  im¬ 
portant  in  recent  years,  and  require  safe  launch 
and  recovery  of  small  boats  and  personnel.  Au¬ 
tonomous  systems  requiring  launch  and  recov¬ 
ery  continue  to  become  more  widely  used  for  a 
variety  of  naval  missions.  Furthermore,  greater 


concern  for  health  and  safety  of  naval  personnel 
has  brought  attention  to  launch  and  recovery, 
which  can  present  increasing  risk  to  crew  as  sea 
state  increases. 

Simulation  is  being  increasingly  used  to  support 
naval  operations  involving  dynamics  of  multi¬ 
ple  bodies.  McTaggart  and  Langlois  [1]  sim¬ 
ulated  replenishment  at  sea  between  a  supply 
ship  and  a  frigate.  Their  simulation  modelled 
the  seaway,  motions  of  both  ships,  and  the  re¬ 
plenishment  gear  used  to  transfer  a  payload 
between  the  ships.  McTaggart,  Roy,  Steinke, 
Nicoll,  and  Perrault  [2]  simulated  launch  and 
recovery  of  a  small  boat  from  a  naval  frigate. 
Their  simulation  modelled  the  seaway,  motions 
of  the  frigate  and  small  boat,  and  the  crane 
used  to  launch  and  recover  the  small  boat. 
Irani,  Kehoe,  Spencer,  Watt,  Gillis,  Carretero 
and  Dubay  [3]  simulated  a  small  unmanned  un¬ 
derwater  vehicle  (UUV)  docking  with  a  subma¬ 
rine. 

Hydrodynamic  interactions  can  significantly  in¬ 
fluence  bodies  that  are  within  proximity  of  each 
other.  Several  authors  have  examined  hydro- 
dynamic  interactions  during  replenishment  at 
sea,  including  McTaggart,  Cumming,  Hsiung, 
and  Li  [4]  and  Andrewartha,  Thomas,  Turner, 
and  Lin  [5].  The  present  paper  examines  hy¬ 
drodynamic  interactions  that  can  occur  dur¬ 
ing  launch  and  recovery  of  a  small  boat  from 
a  frigate. 
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THEORY  FOR  MOTIONS  OF 
A  SINGLE  VESSEL  IN 
WAVES 

Prediction  of  motions  for  a  single  vessel  in 
waves  is  most  commonly  done  using  the  as¬ 
sumption  of  potential  flow,  which  can  give  gen¬ 
erally  accurate  results  with  moderate  comput¬ 
ing  requirements.  The  present  method  is  based 
on  the  work  of  Papanikolaou  and  Schellin  [6], 
which  computes  motions  in  the  frequency  do¬ 
main  using  a  panel  method  and  free  surface 
Green  function  for  zero  forward  speed.  The 
vessel  is  assumed  to  travel  with  quasi-static  for¬ 
ward  speed  and  heading.  This  approach  gives 
generally  good  results  when  the  Froude  num¬ 
ber  Fr  —  U/(g  Lpp )  is  less  than  0.4,  where  U 
is  ship  forward  speed,  g  is  gravitational  accel¬ 
eration,  and  Lpp  is  ship  length  between  perpen¬ 
diculars.  For  sway,  roll,  and  yaw  motions,  lift 
and  viscous  forces  can  be  significant  and  are 
modelled  using  coefficient-based  approaches,  as 
described  by  Schmite  [7]. 

The  motions  of  a  body  with  no  moving  ap¬ 
pendages  can  be  evaluated  in  the  time  domain 
as  follows: 

m  +  [A\)  {77(0}  +  mm)}  +  ic]{vm 

=  {F\t)  +  FD(t)}  (1) 

where  [M]  is  the  ship  mass  matrix,  [A]  is  the 
added  mass  matrix,  is  the  acceleration 

vector,  [B]  is  the  damping  matrix,  {rj(t)}  is 
the  velocity  vector,  [C]  is  the  stiffness  ma¬ 
trix,  is  the  motion  displacement  vector, 

F'(t)}  is  the  incident  wave  force  vector,  and 

F«(t)}  is  the  wave  diffraction  force  vector. 
The  acceleration,  velocity,  and  displacement 
vectors  each  have  six  terms,  representing  surge, 
sway,  heave,  roll,  pitch,  and  yaw  motions.  The 
displacements  rj(t)  in  the  time  domain  can  be 
expressed  in  terms  of  frequency  domain  values 


as  follows: 

rj(t)  =  Real  (rj  e* Ue  (2) 

where  rj  is  complex  motion  amplitude  in  the  fre¬ 
quency  domain  and  uoe  is  wave  encounter  fre¬ 
quency.  The  equation  of  motions  in  the  fre¬ 
quency  domain  can  thus  be  written  as: 

([M]  +  [A])  +  iue[B]  +  [c*])  {rj} 

=  { F 1  +  Fd}  (3) 

Note  that  the  wave  excitation  forces 
F1  +  Fd^  in  the  frequency  domain  are 
complex  quantities,  containing  phase  informa¬ 
tion. 

A  real  ship  will  include  appendages  that  will  in¬ 
troduce  additional  degrees  of  freedom  to  Equa¬ 
tions  (1)  and  (3).  These  degrees  of  freedom 
can  include  rudder  rotations,  stabilizer  fin  rota¬ 
tions,  azimuthing  propeller  rotations,  and  fluid 
motions  within  stabilizer  tanks.  Lloyd  [8]  de¬ 
scribes  incorporation  of  appendages  into  ship 
equations  of  motion. 

THEORY  FOR  MOTIONS  OF 
TWO  VESSELS  IN  WAVES 

The  approach  used  for  computing  the  motions 
of  a  vessel  in  waves  has  been  extended  to  con¬ 
sider  the  presence  of  both  vessels.  Van  Oort- 
merssen  [9]  was  among  the  first  to  compute 
interaction  effects  between  two  bodies  using  a 
panel  method.  Subsequent  related  work  has 
been  performed  by  others,  including  Li,  McTag- 
gart,  and  Hsiung  [10].  When  solving  the  radia¬ 
tion  problem  for  two  bodies,  it  is  necessary  to 
consider  flow  arising  from  motion  of  each  body 
while  the  other  body  is  present  and  stationary. 
The  positions  of  the  ships  relative  to  each  other 
are  assumed  to  be  quasi-steady.  Solution  of  the 
wave  diffraction  problem  considers  both  bod¬ 
ies  as  stationary,  similar  to  the  wave  diffrac¬ 
tion  problem  for  a  single  body.  For  the  case  of 
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neither  vessel  having  moving  appendages,  the 
equations  of  motion  in  the  frequency  domain 
for  two  ships  A  and  B  are  written  as: 
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where  superscripts  A  and  B  refer  to  the  two  in¬ 
dividual  ships.  Note  that  the  added  mass  and 
damping  terms  for  the  the  pair  of  ships  include 
cross-ship  terms  because  motions  of  one  ship 
will  induce  radiation  forces  on  the  other  ship. 
Degrees  of  freedom  for  appendages  such  as  rud¬ 
ders,  stabilizer  fins,  azimuthing  propellers,  and 
U-tube  tanks  can  be  added  to  the  above  equa¬ 
tion  in  a  manner  similar  to  the  single  ship  case 
of  Equation  (3). 


NUMERICAL 

IMPLEMENTATION 

The  ShipMo3D  ship  motion  suite  [11]  has  been 
extended  to  model  seakeeping  in  waves  for  a 
pair  of  vessels.  ShipMo3D  is  written  mostly  in 
C#  [12],  with  some  portions  of  the  code  writ¬ 
ten  in  C++  [13]  to  give  faster  performance  for 
numerical  computations.  ShipMo3D  uses  the 
Math.NET  numerical  library  [14],  which  pro¬ 
vides  efficient  performance  for  tasks  such  as 
solving  systems  of  linear  equations. 

To  eliminate  the  occurrence  of  irregular  fre¬ 
quencies  with  erratic  results  for  predicted  hy¬ 


drodynamic  forces  [15],  panelled  lids  have  been 
implemented  at  the  waterplane  of  each  ship, 
with  a  specified  boundary  condition  of  zero  ver¬ 
tical  velocity.  This  approach  has  been  imple¬ 
mented  for  both  single  ship  and  dual  ship  cases, 
and  has  been  shown  to  be  very  successful  for 
eliminating  irregular  frequencies. 

Computational  requirements  for  ship  interac¬ 
tion  cases  are  typically  much  greater  than  sin¬ 
gle  ship  cases  when  solving  for  radiation  and 
diffraction  forces.  While  a  single  ship  case  nor¬ 
mally  assumes  lateral  symmetry  when  solving 
radiation  and  diffraction  forces,  this  symmetry 
is  typically  lost  when  considering  a  pair  of  ships. 
Consequently,  the  number  of  unknowns  to  be 
solved  increases  by  a  factor  of  approximately  4, 
and  the  number  of  influence  terms  for  solving 
flow  parameters  increases  by  a  factor  of  approx¬ 
imately  16.  In  consideration  of  these  large  com¬ 
putational  requirements,  care  has  been  given  to 
using  memory  efficiently,  including  minimizing 
the  number  of  times  that  large  arrays  are  ini¬ 
tialized  in  memory. 

The  ShipMo3D  application  RadDifPair  com¬ 
putes  a  radiation  and  diffraction  force  database, 
which  is  intended  to  cover  all  combinations  of 
encounter  frequency,  ship  speed,  relative  sea  di¬ 
rection,  and  wave  frequency  that  a  pair  of  ships 
might  encounter.  Typical  values  for  a  pair  of 
ships  with  lengths  between  100  and  200  m  are: 

1.  encounter  frequencies  uje  —  0.05,  0.10,  0.15, 
. . .,  5.0  rad/s  (100  values), 

2.  ship  speeds  U  —  0,  2,  4,  . . .,  30  knots,  (16 
values) , 

3.  relative  sea  directions  f3s  —  0,  15,  30,  . . ., 
360  degrees,  (37  values), 

4.  incident  wave  frequencies  uoi  —  0.10,  0.15, 
0.20,  . . .,  2.0  rad/s. 

Computation  of  added  mass,  damping,  and 
diffraction  force  terms  for  the  above  conditions 


3 


typically  requires  24  hours  and  6  GB  of  mem¬ 
ory  on  a  desktop  workstation  with  8  dual-core 
processors. 

VERIFICATION  AND 
VALIDATION 

Numerical  verification  tests  were  conducted  to 
ensure  that  the  ship  interaction  software  was 
performing  as  expected.  Motion  predictions 
were  conducted  for  two  ships  far  away  from 
each  other,  with  predicted  motions  showing 
no  significant  influence  from  interaction  effects. 
Motion  predictions  were  also  conducted  for  a 
pair  of  identical  ships  in  close  proximity  and 
aligned  longitudinally.  The  surge,  heave,  and 
pitch  motions  were  identical  for  the  two  ships. 
The  sway,  roll,  and  yaw  motions  were  identical 
in  magnitude  but  opposite  in  phase  between  the 
two  ships. 

Validation  of  motion  predictions  was  conducted 
using  the  experiments  of  McTaggart,  Cum- 
ming,  Hsiung,  and  Li  [4].  During  the  exper¬ 
iments,  the  ship  models  were  constrained  in 
surge,  sway,  and  yaw.  These  constraints  were 
included  in  the  numerical  simulations  by  apply¬ 
ing  high  external  stiffness  to  the  ship  models. 
The  numerical  predictions  very  good  agreement 
with  experiments. 

The  ShipMo3D  ship  motion  software  has  been 
validated  for  a  variety  of  vessels;  however,  it 
hasn’t  yet  been  validated  for  a  RHIB. 

EXAMPLE  MOTIONS  OF  A 
SMALL  BOAT  NEAR  A 
FRIGATE 

Motion  computations  have  been  performed  for 
a  rigid-hull  inflatable  boat  (RHIB)  operating 
in  close  proximity  to  a  naval  frigate,  as  would 


occur  during  launch  and  recovery.  The  RHIB 
is  based  on  a  Zodiac  H935,  with  properties 
given  in  Table  1  when  loaded  with  12  person¬ 
nel.  Figure  1  shows  the  hull  surface  model  for 
the  RHIB,  which  formed  the  basis  for  panelling 
the  hull  for  radiation  and  diffraction  compu¬ 
tations.  Figure  2  shows  the  RHIB  model  for 
motion  simulation,  including  the  wet  panelled 
hull  in  yellow,  the  dry  panelled  hull  in  green, 
and  azimuthing  propellers  from  outboard  en¬ 
gines  in  blue. 

The  generic  frigate  is  a  nominal  ship  that  has 
never  been  built,  and  has  properties  given  in 
Table  2.  Figure  3  shows  the  hull  surface  model 
for  the  generic  frigate,  which  formed  the  basis 
for  panelling  the  hull  for  radiation  and  diffrac¬ 
tion  computations.  Figure  4  shows  the  generic 
frigate  model  for  motion  simulation,  including 
the  wet  panelled  hull  in  yellow,  the  dry  pan¬ 
elled  hull  in  green,  the  propellers  in  blue,  and 
the  bilge  keels,  skegs,  propeller  shaft  brackets, 
and  rudder  in  red. 

Single  vessel  motion  computations  were  per¬ 
formed  for  the  RHIB  and  generic  frigate.  Ship 
interaction  computations  were  performed  for 
the  RHIB  along  the  starboard  side  of  the 
frigate,  with  a  1.5  m  gap  between  the  2  ves¬ 
sels.  All  computations  were  performed  for  the 
upper  limit  of  Sea  State  3,  with  a  significant 
wave  height  Hs  of  1.25  m  and  a  peak  wave 
period  Tp  of  7.5  s.  This  sea  state  is  consid¬ 
ered  representative  of  the  upper  limit  at  which 
launch  and  recovery  operations  might  be  safely 
performed.  Ship  speeds  of  0  to  12  knots  were 
evaluated,  encompassing  speeds  that  would  be 
realistically  used  for  launch  and  recovery  oper¬ 
ations.  For  a  ship  speed  of  12  knots,  the  as¬ 
sociated  forward  speed  Froude  number  for  the 
RHIB  is  0.65,  which  likely  is  sufficiently  high 
to  introduce  challenges  for  the  numerical  pre¬ 
dictions. 

Among  the  various  vessel  motion  parameters 
that  will  influence  launch  and  recovery,  the  fol- 
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Length 

9.35  m 

Beam 

3.00  m 

Displacement 

7242  kg 

Number  of  outboard  engines 

2 

Number  of  personnel  onboard 

12 

Natural  roll  period 

2.4  s 

Length  between  perp. 
Beam 

Displacement 
Number  of  rudders 
Number  of  propellers 
Natural  roll  period 


120  m 
14.1  m 
3713  tonnes 
1 
2 

8.9  s 


Table  1:  Properties  of  Rigid-Hull  Inflatable  Boat 


Table  2:  Properties  of  Generic  Frigate 


Figure  1:  Model  of  Hull  Surfaces  for  Rigid-Hull 

Inflatable  Boat  Figure  3:  Model  of  Hull  Surfaces  for  Generic 

Frigate 


Figure  2:  Ship  Motion  Model  for  Rigid-Hull  In-  Figure  4:  ship  Motion  Model  for  Generic  Frigate 
flatable  Boat 
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lowing  are  likely  among  the  most  important: 

•  vertical  acceleration  of  generic  frigate 
crane  tip, 

•  relative  vertical  velocity  between  generic 
frigate  crane  tip  and  ocean  surface, 

•  vertical  acceleration  of  RHIB, 

•  relative  vertical  velocity  between  generic 
frigate  crane  tip  and  RHIB  attachment 
point, 

•  relative  lateral  velocity  between  generic 
frigate  crane  tip  and  RHIB  attachment 
point, 

•  relative  vertical  displacement  between 
generic  frigate  crane  tip  and  RHIB  attach¬ 
ment  point, 

•  relative  lateral  displacement  between 
generic  frigate  crane  tip  and  RHIB 
attachment  point. 

Vertical  acceleration  of  the  crane  tip  will  in¬ 
fluence  loading  on  the  crane,  particularly  when 
the  RHIB  is  emerged  from  the  water  and  the 
crane  is  supporting  its  full  weight.  The  rel¬ 
ative  vertical  velocity  between  the  crane  tip 
and  the  ocean  surface  will  provide  a  good  in- 
dictor  of  relative  motion  between  the  crane  tip 
and  RHIB  because  the  small  RHIB  will  tend 
to  follow  the  ocean  surface.  Vertical  acceler¬ 
ation  of  the  RHIB  will  place  demands  on  the 
RHIB  crew,  including  coupling  and  decoupling 
between  the  crane  hook  and  RHIB  attachment 
point.  Relative  vertical  velocity,  relative  lat¬ 
eral  velocity,  and  relative  vertical  displacement 
will  contribute  to  the  general  complexity  of  the 
operation.  Relative  lateral  displacement  could 
cause  collisions  between  the  RHIB  and  generic 
frigate. 

Figures  5  to  10  show  polar  plots  of  some  re¬ 
sponses  influencing  launch  and  recovery.  For 


each  plot,  speed  is  given  on  the  radial  axis  and 
relative  sea  direction  is  given  on  the  angular 
axis.  The  color  scale  for  each  plot  is  based  on 
the  range  of  numerically  predicted  values  rather 
than  on  any  operational  limits  that  might  exist. 
For  the  current  configuration  with  the  crane  on 
the  starboard  side  of  the  ship,  a  typical  launch 
and  recovery  scenario  would  have  the  ship  and 
RHIB  travelling  at  10  knots  with  the  nominal 
sea  direction  coming  from  the  port  bow  quar¬ 
ter. 

Figure  5  gives  vertical  acceleration  at  the  crane 
tip  on  the  frigate.  The  crane  is  located  amid¬ 
ships  on  the  frigate;  thus,  vertical  accelerations 
are  greater  in  the  vicinity  of  beam  seas,  which 
will  include  roll  and  associated  vertical  motions 
on  the  crane.  Figure  6  shows  the  vertical  ve¬ 
locity  of  the  crane  tip  relative  to  the  ocean  sur¬ 
face.  Note  that  wave  radiation  and  diffraction 
have  need  been  included  when  evaluating  the 
ocean  surface  elevation.  Figures  7  and  8  show 
heave  acceleration  of  the  RHIB  on  its  own  and 
when  near  the  frigate.  Hydrodynamic  interac¬ 
tion  effects  somewhat  influence  heave  acceler¬ 
ations  for  the  RHIB.  Figures  9  and  10  show 
roll  of  the  RHIB  on  its  own  and  when  near 
the  frigate.  Hydrodynamic  interaction  effects 
have  a  dramatic  influence  on  roll  of  the  RHIB, 
especially  in  the  vicinity  of  head  seas.  Note 
that  the  experiments  by  McTaggart,  Cumming, 
Hsiung,  and  Li  [4]  revealed  a  similar  physical 
phenomenon,  in  which  the  presence  of  a  sup¬ 
ply  ship  induced  significant  roll  motions  on  a 
nearby  frigate  in  head  seas.  The  high  roll  mo¬ 
tions  of  the  RHIB  at  higher  speeds  in  the  vicin¬ 
ity  of  head  seas  are  also  due  to  wave  encounter 
periods  being  near  the  RHIB  natural  roll  pe¬ 
riod  of  2.4  s.  Due  to  the  large  magnitude  of  the 
roll  motions  in  Figure  10,  they  would  be  signifi¬ 
cantly  influenced  by  nonlinear  effects;  however, 
these  nonlinear  effects  haven’t  been  modelled 
in  the  present  frequency  domain  computations. 
It  is  recommended  that  time  domain  methods 
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Figure  5:  RMS  Vertical  Acceleration  of  Crane  Tip 
in  Upper  Sea  State  3 


Following  seas 


Figure  7:  RHIB  Heave  Acceleration  in  Upper  Sea 
State  3  without  Interaction 


Following  seas 


Following  seas 


Figure  6:  RMS  Relative  Vertical  Velocity  between 
Crane  Tip  and  Ocean  Surface  in  Upper  Sea  State 
3 


Figure  8:  RHIB  Heave  Acceleration  in  Upper  Sea 
State  3,  1.5  m  Gap  between  Generic  Frigate  and 


RHIB 
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be  investigated  for  modelling  these  large  ampli¬ 
tude  motions. 
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Figure  9:  RHIB  Roll  in  Upper  Sea  State  3  without 
Interaction 


Head  seas 


Following  seas 


Figures  5  to  10  suggest  that  following  seas  could 
be  an  excellent  heading  for  launch  and  recov¬ 
ery;  however,  steering  can  be  challenging  in  fol¬ 
lowing  seas.  Furthermore,  no  results  have  been 
presented  for  other  quantities  which  might  sig¬ 
nificantly  influence  launch  and  recovery  oper¬ 
ations.  Greater  understanding  of  launch  and 
recovery  operational  limits  would  assist  in  de¬ 
termining  which  parameters  should  be  investi¬ 
gated  in  greater  detail. 


CONCLUSION 

Frequency  domain  predictions  of  vessel  motions 
in  waves  have  been  extended  to  predict  mo¬ 
tions  for  two  vessels,  including  hydrodynamic 
interaction  effects.  Example  computations  for 
launch  and  recovery  of  a  RHIB  from  a  frigate 
indicate  that  the  presence  of  the  frigate  can  in¬ 
duce  large  motions  on  the  RHIB.  The  high¬ 
est  RHIB  motions  occur  at  higher  speeds  in 
the  vicinity  of  head  seas,  for  which  wave  en¬ 
counter  periods  are  near  the  RHIB  natural  roll 
period.  For  the  large  amplitude  RHIB  motions 
observed  in  the  example  computations,  it  is  rec¬ 
ommended  that  time  domain  methods  be  inves¬ 
tigated  for  modelling  associated  nonlinear  ef¬ 
fects.  Validation  of  the  RHIB  motion  predic¬ 
tions  is  recommended  due  to  the  high  Froude 
numbers  for  the  RHIB. 
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